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It is known that the transfer of spin angular momenta between current carriers and local moments 
occurs near the interface of magnetic layers when their moments are non-collinear. However, to 
determine the magnitude of the transfer, one should calculate the spin transport properties far 
beyond the interface regions. Based on the spin diffusion equation, we present a self-consistent 
approach to evaluate the spin torque for a number of layered structures. One of the salient features 
is that the longitudinal and transverse components of spin accumulations are inter-twined from 
one layer to the next, and thus, the spin torque could be significantly amplified with respect to 
treatments which concentrate solely on the transport at the interface due to the presence of the 
much longer longitudinal spin diffusion length. We conclude that bare spin currents do not properly 
estimate the spin angular momentum transferred between to the magnetic background; the spin 
transfer that occurs at interfaces should be self-consistently determined by embedding it in our 
globally diffuse transport calculations. 

PACS numbers: 72.25.-b, 72.15.Gd, 73.23.-b 

I. INTRODUCTION 

The concept of using a spin polarized, current to switch the orientation of a magnetic layer was developed by Slon- 
czewski and Bergerd, and has been followed up by several otherso Recently, we have proposed a way to understand 
this spin transfer torque by adapting the model we used to understand magnetoresistance for currents perpendicular 
to the plane of the layer (CPP)B. Namely, two phenomena, CPP magnetoresistance (MR) and spin torque, originate 
from the spin accumulation. The former is primarily associated with the longitudinal spin accumulation and the 
latter is governed by a transverse effect. The distinguishing feature between previous treatments and the one we 
recently outlined lies in our focus on the spin transport for the entire CPP structure rather than for the interface 
region alone. The specular scattering of the current at interfaces between magnetic and nonmagnetic layers that is 
attendant to ballistic transmission can create spin torqueM. Here we start at the opposite extreme and consider the 
spin torque due to the bulk of the magnetic layers and the diffuse scattering at interfaces; as is the case for giant 
magnetoresistance (GMR) reality is probably a mixture of these extreme positions. 

To understand the significance of our approach, we recall the physics of CPP-MR. The resistance of the entire 
CPP structure comes from various sources of scattering. If one only considers scattering at an interface, one may 
concentrate on the calculation of the transmission coefficient for the interface. However, ballistic transmission across 
an interface is not the only physics of the transport in magnetic multilayers. The leads, as well as impurity scattering 
in the layers, have to be included in the calculation of the resistance of the entire CPP structure. Therefore, one 
should embed the ballistic interface scattering in the framework of the diffusive scattering from the bulk of the layers, 
as well as the diffuse interface scattering, for a better approach to describing the transport. This is very much 
the spirit of CPP transport: the spin transport is described by macroscopic spin diffusion while the detail of the 
interface scattering is treated as a boundary ppndition. The crossover from ballistic interface scattering to diffusive 
scattering has been recently studied in detailD. By analogy, we argue that the spin torque should be calculated by 
solving the transport equation for the entire structure. The interface ballistic spin transfer, which is the center of 
the previous discussions, may be physically important; however it should be embedded in a larger picture, i.e., the 
ballistic transport across an interface should be connected to the diffusive transport outside the interfacial region. We 
will show in this paper that our semi-classical formalism supplies a natural framework to incorporate this. Indeed 
Stiles and Zangwill have done just this; however we make a key assumption, that is different from theirs. That is, 
we assume that a component of the spin current transverse to the magnetization exists in the magnetic layersE We 
should emphasize that our transport calculation does not conflict with the physics of ballistic transport at interfaces. 
In fact, we will see below, the ballistic component can be accommodated within our formalism. 
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FIG. 1: Multilayered pillar-like structure used for current induced reversal of a magnetic layer. FM2 is a thick ferromagnetic 

cos 9e z — sin 6e y , Sp is a thin nonmagnetic spacer, FM1 



layer with the thickness exceeding Af^; and local magnetization Mj 



is a thin ferromagnetic layer with the thickness tp and local magnetization M^ 1 ' = e 2 , and NM is a nonmagnetic back layer 



Here we specify a model system to calculate the spin torque: a magnetic multilayer whose essential elements consist 
of a thick magnetic layer, whose primary role is to polarize the current, a thin magnetic layer that is to be switched, 
a nonmagnetic spacer layer so that there is no interlayer exchange coupling between the thick and thin layers, and 
a nonmagnetic layer or lead on back of the thin magnetic layer; see Fig. ^. As we show in this paper the angular 
momentum transferred to a thin layer far exceeds the transverse component (to the orientation of the magnetization 
of the thin layer) of the bare portion of the incoming spin polarized current, i.e., that part proportional to the electric 
field; see Eq. (|) below. It is a direct consequence of the spin accumulation coming from the two primary layers, the 
thick magnetic and nonmagnetic back layers, that produce this buildup. The role of this accumulation in the spin 
current is given in Eq. (7) of Ref. |5|; it is a consequence of considering the transport in the multilayer as a diffusive 
process, and is in keeping with our previous treatments of transport in magnetic multilayers with the proviso that 
one has to include the exchange interaction between the accumulation and the magnetic background, a.k.a. the "sd" 
interaction, to obtain the angular momentum transferred. It is this interaction which produces the spin transfer 
between the current and the magnetic background. Among other things the parameters entering our theory are 
determined from CPP transport measurements, except for the exchange interaction between the itinerant electrons 
and the magnetic background. 

In this paper we first review the formalism presented in Ref. || for calculating the torque and effective field acting on 
a magnetic layer. In particular we define the boundary conditions between the layers of the multilayer and point out 
what are the sources for the longitudinal and transverse spin accumulations when the magnetization of the layers are 
noncollinear. In Sec. Ill we present our results for the spin accumulation, spin currents and the torque and effective 
field acting on the thin layer of the multilayer depicted in Fig. |l|. The bulk of our results are obtained numerically, 
however in certain limits we are able to give analytic expressions, e.g., for the amplification of the torque and effective 
field acting on a layer due to the accumulation. In Sec. IV we present an analytic expression for the amplification, and 



we indicate how recent "ab-initio" determinations of the change of spin currents across a ferromagnetic-nonmagnetic 
(FM /NM) interface can be incorporated at the interfaces between layers in our diffusive treatment of the transport in 
a multilayered structure. The impact of transverse accumulation on GMR is given in Sec. [v[ in particular in creating 
deviations from the simple linear dependence of the resistance with the cosine of the angle between the thick and thin 
layers. 



II. REVIEW OF FORMALISM 



By starting with the linear response of the current to the electric field for diffusive transport! 



j(x) = CE(x) - D^-, 



(1) 
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we find the spin or magnetization current is 

j m = ReTr(crj) 
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where both current j(x) and field E{x) are directed along the growth direction x of the multilayeri. Here C denotes a 
conductivity and D a diffusion constant; the two are related by the Einstein relation C = e 2 N(ep)D for a degenerate 
metal, where N(ep) is the density of states at the Fermi level. The first term on the right hand side of Eq. (||) is 
the bare contribution to the spin polarized current from the electric field, the third term is the contribution from 
the spin accumulation attendant to the current across a magnetically inhomogeneous structure; although it can be 
subsumed as an interface contribution, as we include the diffuse interface scattering below, we will neglect the second 
term coming from charge accumulation. As we will show in magnetic multilayers the contribution of the third term 
to the spin current can dominate over the first. 

The equation of motion for the spin accumulation is 



9m d\ m , T ... „ „ m 

1u + ~a + (-Vfc)m x M d = , 
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(3) 



where t s j is the spin-flip relaxation time of the conduction electron. While the last term is diffusive the spin motion 
induced by the exchange interaction J is not at all diffusive in spin space; it describes a deterministic or ballistic 
rotation of the accumulation that itself is generated by diffusive spin-flip processes. The diffusion equation for the 
macroscopic variables m and j m can be derived from the Boltzmann equation for the spin distribution function in the 
limit where the length scale dj = vphJJ > X m f p , i.e., when the distance an electron moves while its spin rotates by 
2-k is greater than the mean free pathta. In the opposite limit we are not able to derive such a simple relation between 
these variables from the Boltzmann equation; nonetheless we will calculate the spin torque by using this expression 
and later evaluate the significance of the results. We note that in this treatment we consider the effect of the exchange 
interaction between the itinerant electrons and the background magnetization Hint = -Jm • M<j, a.k.a. the "sd" 
interaction, in the equation of motion for the distribution function. To write the spin current Eq. (0) in terms of the 
electric current 



Re(Trj) = 2C E{x) - 2D { 



dn 
dx 



2D 



9m 

dx ' 



we insert this expression in Eq. (H) to eliminate the electric field and charge density, and we obtain 
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where we have dropped an uninteresting term proportional to the derivative of the charge accumulation dn^/dx. Upon 
placing this expression in Eq. (|^) we find the equation of motion for the spin current is 



1 9m 9 2 m • / , <9 2 m\ m m x 1 d 

21* "ft - aF+Pfi™* + ^ + = -™ dx- {fyM (6) 

where /3,j3' are the spin polarization parameters defined by the relations C = /JCoM^, where is the unit vector 
to represent the direction of the local magnetization and D = /3'_D Mj3, and we have defined X s f = ^/2DoT s f and 
Aj = \J2TiDqI J = Xmfpdj/SiT where r s f is the spin flip relaxation time, and J the exchange between the itinerant 
electrons and the magnetic background. The diffusion constant is Dq ~ (l/3)vpT = (l/3)ui?A m / p , where vp is the 
Fermi velocity, r the momentum relaxation time, and X m f p the mean free path; X m f p = Vp r. As we indicate later, 
if one interprets A m / P as that associated with the diffusive scattering at interfaces one arrives at a much smaller 
estimates for Aj than when one uses the mean free path arriving from scattering in the bulk of the layers. 

The term on the right hand side of the time dependent-diffusion equation Eq. (^) for the spin accumulation is the 
source term; it is this term that drives the accumulationlMx Here we will look for the steady state solutions so that 
the first term on the left hand side is zero, and the electric current j e is constant throughout the multilayer. In 
Appendix [A] we discuss the source term; here we point out that this term guarantees the continuity of the spin current 
j m across the interfaces, provided the accumulation is continuous. This can be immediately verified by integrating 
Eq. (Q) across an interface and using the definition of spin current Eq. (Q), or by integrating Eq. (^|). 

We will proceed along the lines of the conventional treatment for current perpendicular to the plane of the layers 
(CPP) in magnetic multilayers and focus on the discontinuous variation of the background magnetization between 
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the layers. In this treatment of CPP we assume thai, the magnetization is uniform throughout a layer so that the 
source term is confined to interfaces between layersEJ,E£l; in this case one can take into account the source terms by 
appropriate boundarjL.conditions; this is the procedure usually followed when calculating the spin accumulation in 
magnetic multilayersEj. In this case we set the source term in Eq. (^) to zero, separate the spin accumulation into 
longitudinal (parallel to the local moment) and transverse (perpendicular to the local moment) modes, and look for 
the steady state solutions. We stress that the terms longitudinal and transverse are relative to the magnetization in 
the individual layers, i.e., they are locally defined and have no global meaning throughout a multilayer. Equation (^ 
) can now be written as 

<9 2 mii mil , 
- = 0, (7) 



where X sd i = y/1 - (3[3'X s f, and 



9x2 X sdl 



d 2 mj^ m | m x M d 



(8) 



The longitudinal accumulation rxi|| decays at the length scale of the spin diffusion length X sd i while the transverse spin 
accumulation mj_ decays as Aj. It is important to differentiate between X sd i which arises from spin flip processes, 
and A j which represents the decay of transverse spin currents due to ordinary spin dependent, \ m f p , but non-spin- flip 
scattering. For cobalt we estimate in Appendix |A| the transverse spin accumulation has a much shorter length scale 
compared to the longitudinal one; for permalloy Aj is comparable to the spin diffusion length. We limit our present 
study to A,/ < Xsdi ■ 

The boundary conditions are that the spin accumulation m and current j m are continuous across at interfaces as 
long as there is no specular or diffusive interface scattering; note that it is unnecessary to invoke the delta function 
source term at the boundary of a layer (see the Appendix |A|) which specifies the discontinuity in dm/dx. As noted 
above the source of this discontinuity in the accumulation is to guarantee the continuity of the spin current (provided 
no torque is created at interface, see below), therefore we can simply invoke it; see Appendix |^. While we could in 
principle .consider the effect of specular scattering at interfaces, as we have for other problems associated with CPP 
transportE3, this involves considering each k vector separately which for a noncollinear multilayer is a cumbersome 
problem we will not address. Othcjt&.have considered specular scattering at interfaces and have shown it gives rise 
to spin torque at these boundariesElcl. The presence of diffuse scattering at interfaces due to both roughness and 
interdiffusion in noncollinear structures is treated in Appendix |b[ 

Once having the spin accumulation we take a look at its influence on the background magnetization. The equation 
of motion for the local magnetization is 

— - =-7oM d x H e + Jm )+aM d x — -2, 9 
at dt 

where 70 is the gyromagnetic ratio, H e is the magnetic field including the contributions from the external field, 
anisotropy and magnetostatic field, the additional effective field Jm is due to coupling between the local moments 
(background magnetization) and the spin accumulation, and the last term is the Gilbert damping term. As seen 
from Eq. (|), the longitudinal spin accumulation has no effect on the local moment; we can re- write Eq. (^) in terms 
of the transverse spin accumulation only by replacing m by mi. As shown in Ref. 3 the two components of the 
accumulation in the plane transverse to the magnetization M^ 1 ' of the layer for which we are calculating the effect 
due to the spin current are 

Jm ± = aM|j 2) x + bM^ (10) 

where is the magnetization of the other layer; in the case we discuss in this paper (see Fig. [I]) M^ 1 ' refers to 

the thin layer which is being switched and to the thick layer which polarizes the current and which is pinned 

so that is does not rotate. By placing this form of the accumulation in the equation of motion for the background 
magnetization Eq. (JsJ) we find that the transverse spin accumulation produces two effects simultaneously: the term 
bM^x m|j 2) is the torque due to an "effective field" bM^ , and the other is aM^ x (M^ x M^ 1 ') which is called 
the "spin torque" predicted by Slonczewski and BergerB. The first term produces a precessional motion about M^ 1 ''; 
in this sense it acts as if the spin current creates a magnetic field on mJj 1 '. The second term acts so as to increase or 
decrease the angle between Mj and M^; also, it acts so as to assist or oppose the damping term in Eq. (||). We 
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stress that it is the "sd" exchange interaction between the spin accumulation attendant to CPP and the background 
magnetization that are the origins of the spin torque a and effective field b . 

Another way of determining the torque transmitted by the current to the background comes from recognizing that 
angular momentum is conserved so that m + = const, and 

t ~ dm/dt = -dMd/dt. (11) 

By following this alternate path, which is indeed the path taken by most who have worked on this problem, we write 

r ~ dm/dt + djrn/dx, (12) 

which says that the torque transmitted by a steady state current is given by the gradient of the spin current. When 
one integrates this over a layer, or even across an interface, which absorbs the momentum we find 



At = 



If 



(dj m /dx)dx =j m (t F ) - j m (0). 



(13) 



But from Eqs. (§), (|) and (§) 



(dim/dx) ~ -2D (d 2 m/dx 2 ) ~ -2D ^ 2 d . (14) 

From this form it follows that the x component of the torque comes from the y component of the accumulation as 
Md is defined as the z direction in the layer which is receiving the spin angular momentum. So this is contrary what 
one may think from hastily looking at Eq. (JbJ) where one sees that the x component of the spin current is related 
to the x component of the gradient of the spin accumulation; however for the torque it is-ithje gradient of the spin 
current that enters. One can also apply Eq. ([l3]) across an interface as ab-initio calculations! 2 ) □ have shown that spin 
dependent specular reflections indeed induce a torque at interfaces. In this case Eq. ([l4]) should be interpreted as the 
discontinuity of the spin current at an interface; this is concomitant to either a rotation of the accumulation to across 
the interface, a discontinuity in to, or a discontinuity in the background magnetization as we discuss in Appendix [a|. 

In the next section we present our results for the spin accumulation and current, and the ensuing spin torque and 
field acting on the thin magnetic layer for the multilayer depicted in Fig. |l| 



III. RESULTS 

The essential elements of the multilayered pillar- like structure used for current induced reversal of a magnetic laye 
are shown in Fig. The nonmagnetic lead in back of the thick magnetic layer is not necessary for our discussion 
as we have taken the thickness of the ferromagnetic layer to the left to exceed X s di', what matters is that the spin 
polarization of the current is primarily dictated by the thick magnetic layer. Depending on the thickness of the thin 
magnetic layer it can have an effect on the longitudinal accumulation, however for the thicknesses tp of interest its 
role is minor. We have to solve for the spin accumulation m and spin current j m in four layers with 3 interfaces: 1) 
at the thick magnetic layer and the nonmagnetic spacer layer, 2) between the spacer and the thin magnetic layer, 
and 3) between the thin layer and the nonmagnetic back layer. At each interface there are 3 components each for m 
and j m to match; in all 18 parameters. One approximation, which is valid for all cases of interest, is to consider the 
nonmagnetic spacer layer thickness small compared to the spin diffusion length in the nonmagnetic spacer, i.e., tjv 
■C X^ d i ~ 600 nm. Then both m and j m are constant across tjy in which case we can focus on three layers and only 
12 parameters. In Appendices [b| and |c] we derive the boundary conditions on the accumulation and current; with 
these we can determine the spin current across the entire structure, and consequently the spin torque and effective 
field acting on the thin magnetic layer. Without further simplifications we are unable to give analytic expressions for 
the accumulation and current across the multilayer, and present our numerical results for these quantities as well as 
the torque and field they create. In the following section we derive an analytic expression in a limiting case. 

In Fig. H we show the total (not per unit length as in Ref. ||) spin torque and effective field, see Eq. (|l0|), as a 
function of the thickness of the thin magnetic layer tp which is being switched; for these figures we have taken A,/ = 4 
nm, which is comparable to \ m fp = 6 nm in the bulk, and \f dl = 60 nm. While this value of Aj may be larger 
than what one should use for say Co the plots clearly indicate the new phenomena that occur around the interfaces. 
In these plots we have considered neither specular nor diffuse scattering at the interfaces; the diffusion constant Dq 
is taken to be 10 -3 m 2 /s in the magnetic layers, and 5T0 -3 m 2 /s in the nonmagnetic layers here, as well as in all 
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FIG. 2: Total torque a sin Otp/Pje ^j^- (a), and total effective field bsaxOtF / Pje-^- (b) acting on the thin ferromagnetic layer 
as a function of its thickness ti? for Aj = 4 nm, Af d! = 60 nm, and zero interface resistance ARi — 0. 

following plots. While the torque rapidly increases for small but finite tp Aj and then gradually levels off, the field 
is largest about tp w 0.5Aj and then decreases towards zero with tp; this can be understood as follows. When the 
thickness of the thin layer tp is much smaller than Aj the spin accumulation in the thin layer is the same as that 

(2) 

of the thick layer at the interface, and its direction is parallel to that of the magnetization of the thick layer Mj 
(remember we are not considering torques created at interfaces); therefore, only the effective field exists, see Eq. (JTc|) . 

(2) — 

As tp increases the spin accumulation in the thin magnetic layer rotates away from and develops a transverse 

component x M^ 1 '', i.e., a spin torque develops. Indeed when tp A,/, the spin accumulation is rotated in the 

thin layer and thus the component of the spin accumulation in the plane of the magnetizations decreases rapidly, i.e., 
the effective field diminishes faster than the torque. As tp increases further there are no additional contributions to 
either the field and torque because they represent effects that are centered at the interface with the spacer layer and 
averaged over the entire thickness of the thin magnetic layer. Although the effective field is negligible compared to 
the torque in the limit of large tp it is noteworthy that at its maximum the field b is at least as large as a. One also 
notes that while the torque and field terms a and b are largest for 9 = 180° and 0° they do not act on the background 
magnetization because sin 9 = 0; here 9 is the angle between the magnetizations of the thick and thin magnetic layers. 
The largest effects are found for 9 ~ 150° - 170°. 

In Fig. ^ we show the accumulation and spin current that produce these torques and fields for a thin layer thickness 
tp = 3 nm which is close to where the torque term starts to saturate. As both m and j m are constant across 
tjy we do not show the nonmagnetic spacer layer in these plots, i.e., we plot m and j m as if tjv = 0. Also we 
take the magnetization in the thin layer M^ 1 -* as our global z axis; the current is along the x axis which is along 
the growth direction of the multilayer, and the y direction is perpendicular to the other two. In these global axes 
longitudinal and transverse in the thin magnetic layer refers to the directions z and x — y; however for the thick layer, 
whose magnetization M^ 2 ' is at an angle 9 relative to M^ 1 ', the global y and z axes do not define what is meant by 
longitudinal and transverse in this layer. Also in the nonmagnetic layers, where there is no equilibrium magnetization, 
we talk only about longitudinal accumulation. Far from the interface x <C — Aj the accumulation and current in the 

(2) 

thick layer are collinear with background magnetization , i.e., referred to its local axes they are longitudinal with 

no transverse components, and the spin current approaches its bare value j m — > /3j e 'M.^ (see Eq. (||)); even though 
one still has a longitudinal spin accumulation in the region — A^ ; <C x <C — Aj its gradient is small compared to that 
of the transverse accumulation which makes large contributions to the spin current in the region x > —Xj. This 
is clear from the plots for m x and j m ,x which goes to zero, while m y — > msin#, j m y — > /3j e sin(9,m z — > mcos(9, and 
jm.z — * Pj e cos9. With this identification it becomes clear for example why for 9 = 90°, j m<y — > 1, in units of (3j e , 
while j m ,z -> 0. 

The results for the spin current in Fig. ^ are interesting; far to the left the current is polarized along as 
one expects in the bulk of a ferromagnetic. Also in the nonmagnetic layer to the right x > tp there is no, or very 
little, spin-cuipent in the x — y directions; the current is polarized along i.e., the region — Aj < x < tp acts as 

a "spin filter"li3 inasmuch as the component of the current transverse to the magnetization of the thin layer which 
is being switched is "absorbed" within this region. That much has been predicted by most treatments of current 
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induced switchingOH. The surprise lies when we look at the enhanced x — y or transverse components of the spin 
current in the region about x = 0. As the "torques" (what we call the spin torque and effective field) transmitted 
to the thin layer tp are just the difference between the transverse components of the spin current at the boundaries 
of the thin magnetic layer, see Eq. (|l3|), we find t x = j m , x (x = tp) — j m ,x(% = 0) ~ b for the effective field , and 
T y = jm,y( x = £f) — jm,y( x = 0) ~ a for the spin torque are both amplified when compared to what finds when 
one neglects spin accumulation. The thick magnetic layer to the left x < is pinned so that the enhanced torques 
acting in the region of the interface do not produce any rotation. The z or longitudinal component of the incoming 
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spin current is not absorbed by the thin magnetic layer as there is no transfer of spin angular momentum along this 
direction (tp -C X s di)', see Fig. ||f. The slight decrease in j m<z is due to the ambient spin flip scattering in the magnetic 
layer which is characterized by A^ ; ~ 60 nm in the plots shown in Fig. ||. The much slower decrease in j m ,z(x > tp) 
comes from the spin flip scattering in the nonmagnetic layer whose A^ ; ~ 600 nm 3> Xf d[ . 

The large enhancement of the transverse spin currents can be understood as follows. Around x — 0, which in our 
picture contains the interfaces between the thick and thin magnetic layers, the source term for the transverse spin 
accumulation is comparable to that for the longitudinal; as mentioned in the previous section the spin accumulation 
makes up for the discontinuity in the "bare" spin current. At the interface between the thick and thin layers this is 
the component of the spin current, coming from the bulk of the thick magnetic layer, that is perpendicular to the 
magnetization in the thin layer. The distance over which it is absorbed is much smaller than that for the longitudinal 
accumulation Xj -C A^ ; . Therefore the gradient of the transverse accumulation about x — is large and as it is 
the gradient that contributes to the spin current, Eq. (^), we find the unanticipated amplification of the transverse 
components of the spin current at this interface (really two interfaces). This amplification is not a maximum about 
9 = 90°, because in addition to the "bare" contribution, there is the component of the spin current that arises from 
the longitudinal accumulation in the thick layer, i.e., parallel to the magnetization Mj 2 \ For 9 — 90° this longitudinal 
accumulation is quite small, see Fig. ^|h, so that there's little amplification at this angle. In the next section we present 
a more quantitative reason for the enhancement. 

The plots in Fig. || were for the case of A^ ; = 60 nm, A,/ = 4 nm, tp ~ where the torque a starts to saturate. To 
determine the roles of the spin diffusion length, A^ ; (while keeping A^ ; = 600 nm), the spin transfer length Aj, and 
the interface resistance due to diffuse scattering at the interfaces ARi (see Appendix |b|), on these plots we have rerun 
our program for: Xf d[ = 30 nm, A,/ = 4 nm, ARj — 0, to show effect of reduced X s di on torques; Xf d[ = 60 nm, Xj = 1 
nm, ARi = 0, to show effect of smaller Xj on torques; and Xf dl = 30 — 60 nm, Xj = 1 — 4 nm, ARi ^ 0, to show 
effect of interface scattering on torques. We account for this scattering in Appendix |c] by introducing thin "interfacial 
regions" in the magnetic layers which have the enhanced scattering found at interfaces, and derive the boundary 
conditions on the accumulation and current in the presence of interface scattering as we let the thickness of the region 
tend to Zjeco. The amount of the scattering ARi and its spin dependence 7 are taken from experimental data on 
CPP-MRO. While the diffuse interface scattering by itself produces sizeable discontinuities in the accumulation it 
does not create much spin torque; this is one difference between our treatment and others. 

In Figs. ^ and ^| we summarize our findings by showing plots of the spin torque a sin 9 and effective field bsin.9 
as functions of the thickness tp of the thin magnetic layer for different combinations of Aj, Xf d[ , ARj, and for two 
different angles 9, 30° and 150°, between the magnetic layers. We find that interface resistance increases these torques 
and causes the spin torque to achieve saturation for smaller tp. By reducing the spin diffusion length to A^ ; = 30nm, 
we find the spin torque and effective field are reduced. When we reduce the spin transfer length to A,/ = 1 nm we 
find the spin torque achieves saturation for smaller tF and the effective field is increased and peaks for lower tp . 

While the number of plots for the accumulation and currents for different par ameters and thicknesses tp are too 
numerous to be shown in this paper, one can view them, in color, at our website http:/ /physics. nyu.edu/~avs203 . 



IV. SPIN TRANSFER LAYERS 



From the results presented in the last section we arrive at a picture o£..khc region in which spin angular momentum 
is transferred that is somewhat different from the conventional pictureoB, i.e., that the spin transfer, as it is called, 
occurs within the few angstroms or monolayers of the thin magnetic layer that is being switched. In the conventional 
picture we would find the transverse spin current in the ferromagnet would be zero near x 0, and the smaller Aj 
the closer to x = 0. In the thick magnetic layer we find the transverse components of the spin current goes to zero as 
expected from Eq. (|J), however, for the thin layer we find the region of spin transfer, defined as the one in which the 
transverse components of the spin current are absorbed in Fig. |3j, is over the entire thin layer when we take tp = 3 
nm ; this corresponds to £ = tp/\/2Xj = 0.5 for the Xj = 4 nm we chose. For thinner layers, e.g., tp — llA, the 
transverse spin current j m ,y is not entirely absorbed by the layer. To better understand this behavior we varied the 
length scale Xj from 4 nm to 1 nm, while keeping the tp — 3 nm the same, i.e., we varied £ between 0.5 and 2, and 
found the transverse current still goes to zero only at x = tp. However, for larger thicknesses of tp such that £ > 5 
and by reducing the resistivity of the nonmagnetic back layer relative to the magnetic layers (typically by a factor 
of 5), we do find the anticipated behavior, i.e., the transverse spin current in the thin magnetic layer goes to zero 
before reaching the interface with the nonmagnetic back layer x — tp. The reason for the different behavior in the 
thin layer arises from its confined geometry, i.e., the reflections from the thin magnetic layer/nonmagnetic back layer 
interface create the patterns observed for the transverse spin currents. 

On the basis of our results we can model the transfer of the spin angular momentum as occurring over a region 
of several lengths of size Aj in the magnetic layers. We arrived at this result by assuming the transport in the spin 
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FIG. 4: Total torque a sin 8t f/P je ^j^- (a,c), and total effective field bsinOtF / Pje acting on the thin ferromagnetic layer as 

a function of its thickness tF for 6 = 30° for different values of Xj and \ f d i in this layer, with and without interface resistance. 
Case \.j = 4 nm, X sl u — 60 nm is represented by the solid line, case Aj = 4 nm, X s< u = 30 is represented by the short-dashed 
line, case Xj = 1 nm, A s( j; = 60 by the dashed-dotted line, and case Xj = 1 nm, X s di = 30 by the long-dashed line. 



transfer region is diffusive, so that Aj = ^/ \ m f p dj /3n. In a ballistic treatment of this region one would posit that 
whatever transverse spin current enters the thin layer is absorbed by it over a characteristic length dk = 2% \k-\ — ki | 1 , 
i.e., one could replace ouil diffusive results with ab-initio calculations of the spin transfer across the region which is 
within dk of the interfacca Having noticed this difference between diffusive and ballistic treatments of transport 
near the interface, it is still necessary to perform a global diffusive calculation of the spin accumulation attendant 
to CPP across the entire pillar-like structures to determine the actual spin currents and concomitantly the angular 
momentum transferred to the background magnetization; j m is not just its bare value (3j e lS/Ld in these multilayered 
structures. Instead, one should self-consistently determine the transverse and longitudinal spin accumulation by taking 
into account the entire structure as we have shown in the previous section; here we present a simplified calculation 
that gives an analytic expression for the amplification. 

Let us consider the limit that the transverse components of the spin current decay on the length scale dk <S tih^n -W- 
Based on the band structure and the quantum mechanical probability spin current dk is only a few angstromscl We 
can incorporate this picture of interfacial thin "spin transfer regions" into our global diffuse picture of the transport 
in the same way as we inserted scattering at interfaces (specular and diffuse) in our calculation of CPP resistancefij. 
That is, we posit that the spin accumulation in the region to the left of thin spacer layer x < —d, where d = dk in 
the ballistic limit and d — Xj in the diffusive limit, is (see Fig.l) 

m(x) = AM 2 e x/x ° dl , (15) 
while for tN+d<x<t^+tp in the thin magnetic layer that is to be switched, 

m(x) = (Be x / X ^' + Ce- a / A '*)Mi, (16) 
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FIG. 5: Total torque a sin 8t f/P je ^j^- (a,c), and total effective field 6 sin 0tp/j3je acting on the thin ferromagnetic layer as 
a function of its thickness tF for 6 — 150° for different values of Aj and \f dl in this layer, with and without interface resistance. 
Case Xj = 4 nm, X s di = 60 nm is represented by the solid line, case Aj = 4 nm, X s di = 30 is represented by the short-dashed 
line, case Xj = 1 nm, X s di = 60 by the dashed-dotted line, and case Xj = 1 nm, X s di = 30 by the long-dashed line. 



and in the nonmagnetic back layer for x > tjy + tp , 

m(i) = &" l/A ^Mi. (17) 

The constants are obtained through the boundary conditions as we now show. These relations explicitly state that 
the spin accumulation is purely longitudinal once the electrons are in the "bulk" of the ferromagnetic layers (thick 
and thin layers) and outside the interfacial spin transfer regions. To determine the spin torque for each individual 
ferromagnetic layer, one needs to find the spin current jo and spin accumulation mo in the nonmagnetic spacer layer. 
By assuming that the spacer layer tjy is small compared to A^ ; , jo and mo are constant across ■ The transverse 
components of jo and mo vary rapidly in the interfacial regions at —d < x < and ijv < x < + d (indeed they may 
be discontinuous), when compared to the longitudinal spin current and accumulation which are continuous across the 
interfaces between the spacer and the ferromagnetic layers, i.e., are constant across the entire region — d < x < t^ + d; 



therefore the latter provide the bridge between the accumulations in the regions x < —d, Eq. (15) and those for 
x > tpf + d, Eq. (|l6|). For example, 

m • M 2 = A, (18) 

and by placing Eq. ([[5]) in Eq. (||) we find 

jo • M 2 = p 3e - 2D a A/Xf dl (19) 
for the interface at x — 0. Similar expressions can be written down at the x = tjv and x — ijv + tF interfaces. 
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As we postulate that the transverse component of the accumulation and current is limited to a spin transfer region 
of size d at the interfaces we have 

*E£ x Ml , 2 =^ x Ml , 2) (20) 
ax a 

and as the spin current is related to the gradient of the accumulation by Eq. (||) we find that jo x Mi,2=±(-Do/d)mo x 
Mi_2 is approximately valid. By using this relation along with the boundary conditions mentioned above, we arrive 
at the transverse current density in the spacer layer 



Jo = ftUsin 2 9 + (2d/X sf ) cos 2 0] _1 Mi x (Mi x M 2 ) (21) 

where we have taken the limit that d <C tp <C \ s f- Clearly, as 9 goes to zero or ir, the magnitude of the spin 
torque is enhanced by a factor of X s f/2d compared to the "bare" transverse current f3j e M.i x (Mj x M 2 ). This huge 
enhancement comes from interplay between longitudinal and transverse accumulation; it is the result of the global 
nature of the spin current even though the transverse component of the spin current and accumulation are absorbed 
within a region d of the interfaces. One should not take the limit d — > because the assumptions we made about 
the spin accumulation, Eqs. ( p"5fjl^ ) and (|2^) break down, i.e., in our calculation the transverse current jg cannot be 
absorbed within d as it tends to zero, and the enhancement does not blow up. 

There is an analogy with how one treats depletion layers in semiconductor p-n junctions; while the transport is 
treated by the diffusion equation the characteristics of the depletion layers themselves are determined from quantum 
mechanics. Similarly while the matching of the Boltzmann distribution functions across interfaces are described by 
quantum mechanics, the overall transport in the magnetic multilayered structure is a problem of diffusive transport. 



V. CORRECTIONS TO CPP RESISTANCE 



The resistance of a magnetic multilayer for CPP has been extensively discussed. At first one limited oneself 
to nominally-. callinear configurations of the magnetic layers, i.e., ferromagnetic and antiferromagnetic alignments 
9 = 0, 180 c O,E3, and noncollinear structures were considered where we took account only of the spin dependent 
scattering through layer dependent self energies, but left out the. effect of the background magnetization on the band 
structure, i.e., we considered conduction by free electron statesEZl. The effect of band structure on CPP resistance has 
been considered by Vedyayevlla. 

When the idea of current induced switching was proposed it was immediately recognized that the transfer of angular 
momentum from the polarized current would have an effect on the voltage drop across the multilayer being studiedta; 
since that time there have been several calculations of the CPP resistance as a function of the angle between magnetic 
layerscr □. Also there has been experimental data on several multilayered -structures that have confirmed that there 
are corrections to the simple cos 2 (#/2) dependence of the CPP resistanceEl However one impediment was that for 
multilayered structures one does not have a good knowledge of the orientation of the magnetization for the individual 
layers, so that one does not have good data on the angular variation of the resistivity. Recently a study of this was 
made on an exchange biased spin valve (ESBV) so as to have a precise determination of the angular dependence of 
the CPP-MREj. The normalized angular dependence of the resistance was defined as 

R(d)-R(0) 

Unarm - _ R{Q) (ZZ) 

and the data was fit to 

l-cos 2 (fl/2) 

«™ - 1 + xcos 2(0/ 2 )- ^ 

Here we present our calculation of the angular dependence of the CPP resistance based on the spin currents we 
find by using the diffusion equation for the spin accumulation in noncollinear structures; see Sec. ||. By treating two 
"thick" magnetic layers, i.e., neglecting reflections from the outer boundaries of the layers, with the magnetizations 
M j 1 ' = cos(#/2)e z + sui(9/2)e y , M?' = cos(#/2)e z — sin(6'/2)e y , where x is the direction of the electric current, with 
a nonmagnetic spacer tjq <C A^ z , and zero interface resistance we obtain the Eq. (|2^) with 

X = j - 1, (24) 
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where 

(1 - gg)Aj 

A =^sr- (25) 

It should be stressed that this expression for x is based oii-Jhe assumption that Aj <C Af d/ . For cobalt, Aj is 
of the order of 3 nm, A^ ; is about 60 nmn3, taking (3 to be 0.919, and calculating 0' « 0.9 by using the densities of 
states for up and down electrons^, we estimate \ to be about 50. We are unable to compare this to data on systems 
containing cobalt; the one system that has been accurately measured has been a series of Py(t)/Cu(20nm)/Py(t) 
ESBV with variable permalloy thicknesses ranging from 6 — 24 nm. The best fit yielded \ — 1-17; however as 
Aj = \f dl for permalloy our expressions for x are not applicable. 

By taking into account the resistance of the interfaces between two FM layers and the normal metal spacer, ARj, 
we obtain 

1 r 
X = A(l+r) ~ 1+ (l + r)(l- 77 ')' (26) 

where r = ARje 2 NQ(ep)(l — 7 ,,2 )/(l — 77'), e is the electron charge, Nq is the density of states at the interface, 7, 7', 
7" are the spin polarization parameters for the conductivity, diffusion constant, and density of states at the interface 
(see Appendix 2). We estimate x to be about 31 for cobalt. 



VI. DISCUSSION OF RESULTS 



The salient conclusion we arrive at is that the bare currents do not correctly estimate the amount of spin angular 
momentum transferred from the polarized current to the background magnetization of the magnetic layers, in the 
layered structures that have been studied to date. It is necessary to do a complete "globally diffusive" transport 
calculation, with the possibility of interfacial ballistic inserts to account for the spin transfer there, in order to 
ascertain the enhancement of this spin transfer by the accumulation attendant to CPP transport. The size of the spin 
transfer region has not been resolved, but we can circumvent this uncertainty by postulating a region d in which a 
transverse component of the_accumulation and current exist, and we can place in this sector either results obtained 
from ab-initio calculation!] or our diffusive spin transfer. Also when we consider diffuse interface scattering the 
mean free path in the region of the spin transfer is considerably smaller, by at least one order of magnitude, than in 
the bulk of the layers. 

The uncertainty in the size of the spin transfer region comes from estimating the magnitude of the "sd" exchange 
interaction J. If one erroneously identifies it with the spin splitting, A ~ 1 eV found from band structure calculations 
which are limited to the diagonal spin components of the exchange-correlation potential, one would indeed find a 
spin transfer region no larger than dk\ however it has been stressed that J ~ kTc ~ 0.1 eV should be identified as 
the "Heisenberg-like" exchange, coupling found in calculations that include the off-diagonal spin components of the 
exchange-correlation potentials. The one case in which one has been able to directly measure J from transmission 
conduction electron spin resonance one found J = 0.106 eV (see Appendix |a|) for permalloyEj; unfortunately no data 
exists for Ccc3. Indeed when we use this value for J and the \ m f P in the bulk of Co we find A,/ ~ 3 nm in which 
case the spin transfer region would be larger than the dk of the order of several angstroms anticipated by othersa □. 
However, when we use the A m / P appropriate to the interface between Co and CuE3, we find Aj ~ 1 nm which is 
comparable to the upper limit estimated from data on Co/Cu pillarsBj. 

Our conclusion about the amplification of the spin torque is independent of the size of the spin transfer region (as 
long as d is large enough so that we can consider the conduction in the semi-classical approximation) because our 
overall calculation of the diffusive transport outside the spin transfer regions remains valid, inasmuch as it is identical 
to the well established theory of Valet and Fert for CPP transport. However the magnitude of the amplification of the 
spin torque does depends on d; see Eq. (|2l|). The size of the spin torque transmitted to the background is primarily 
governed by: the spin dependent transport parameters of the thick magnetic layer which creates the polarized spin 
current, j3, A^, the spin dependent interface scattering parameter 7, and the resistivity and spin diffusion length in 
the normal back layer relative to that of the thick magnetic layer. The characteristics for the relatively thin magnetic 
layers tp used to observe current induced switching do not determine the overall spin accumulation and current in 
the sample; other than sensing the spin polarized current through the "sd" exchange interaction they do not affect 
the size of the spin torque acting on the thin layer. 

In our treatment of current induced switching we considered spin transport across the entire CPP structure rather 
than for the interface region alone, i.e., we have considered the spin torque due to the bulk of the magnetic layers and 
that arising from diffuse scattering at interfaces. We have not considered the problem of matching the distribution 
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functions across adjacent layers in the presence of specular scattering at the interface; this requires a knowledge of the 
band structure in these layers and is outside the scope of our study. The parameters entering our theory arc determined 
from CPP transport measurements, except for J, the "sd" exchange interaction. Previous treatments highlighted the 
spin torque that is attendant to ballistic transmission across an interface between magnetic and nonmagnetic layers; 
as is the case for GMR reality is probably a mixture of these two different positions. 
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APPENDIX A 



We separate the spin accumulation into longitudinal (parallel to the local moment) and transverse (perpendicular 
to the local moment) modes. Equation (^|) can now be written as 

d 2m ll m ll _ 3e dp M , An 



where X s di = yl — pp'\ s f, D[ = ^1 — fifi'Do and p(x) = /JM^, which represents the bare spin polarization of the 
current coming solely from the electric current in the absence of spin accumulation (see Eqs. (^) and ((^)), and 



(9 2 m_L nix m x M a 



dx 2 



\ 2 

A sf 



2D 



dp 

dx 



x M rl 



dp 

dx 



x M d x M„ 



(A2) 



The longitudinal accumulation ni|| decays at the length scale of the spin diffusion length X s di while the transverse 
spin accumulation mj_ decays as Aj. In-a typical transition-metal ferromagnet, e.g. Co, the spin diffusion length A sc y 
has been measured to be about 6Q-JimE3. We estimate Aj by taking the typical diffusion constant of a metal to be 
3 ■ 10~ 3 m 2 /s and J — 0.1 — 0.4 eVc3 so that Xj is about 1.5 nm to 3 nm. Thus, the transverse spin accumulation has 
a much shorter length scale compared to the longitudinal one; it is larger than the mean free path in the interfacial 
region between Co and Cu, ~ 1 nm, and is comparable to X m fp ~ 6 nm in thq-bulk of Co. For permalloy, where we 
can use the value of J = 0.106 eV measured by conduction electron resonances Xj ~ 3 nm which is comparable to 
Am/p and X s di for permalloy. Therefore for multilayers containing Py our treatment is not directly applicable as we 
assume in most of our work Xj <C X sc ii ; we have to go back to the diffusion equation in Sec. [ij and solve the equations 
in this limit. 

It makes a difference whether one treats the magnetization Md(x) as a continuous function or as a finite difference. 
For example, in a domain wall, where one treats the magnetization as a continuously rotating vector, there is no 
longitudinal component of the spin accumulation ni|| coming from the interior of the wall itself because -J^p ~^(M,j) 
is perpendicular, "tangential", to M<j. In this case the transverse component of the source term exists. In our 
treatment of transport across a domain wall we accounted for the continuously rotating magnetization in the wall 
by determining the correction to the electron states induced by the rotations in spin spaceE3; we did not consider 
any spin accumulation. Another treatment of the same problem by SimanekE3 took an approach for the domain wall 
which is more consonant with the equation of motion method we follow in this paper. In that approach the transverse 
spin accumulation due to the continuously rotating magnetization was determined, and he was able to calculate its 
contribution to the domain wall resistance. 

This continuous treatment for domain walls has to be contrasted with the conventional treatment for current 
perpendicular to the plane of the layers (CPP) in magnetic multilayers where one focuses on the discontinuous 
variation of the magnetization between the layers; in this case one indeed does have a longitudinal source term for 
the spin accumulation as we now show. In the usual treatment of CPP we assumq-t|bat the magnetization is uniform 
throughout a layer so that the source term is confined to interfaces between layersciljEJ; in this case one can take into 
account the source terms by appropriate boundary conditions; this is the procedure usually followed when calculating 
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the spin accumulation in magnetic multilayersS By discretizing the source terms in Eqs. (Al) and (A2) we find for 
the longitudinal accumulation 



2Dn 



Miil-Mi-MMxj), 



j=i±i 

while for the transverse accumulation the source term is 

Pi 



2D, 



^ x {Mi x Mj)5( Xj ), 



(A3) 



(A4) 



j=i±i 



where the layer we are considering is labelled i while the interfaces with the adjacent layers j = i ± 1 are at Xj . For 
collinear structures we see that the transverse term is zero; the longitudinal source term at the interface between two 
identical magnetic layers is zero if they are parallel, and 



L>0 



(A5) 



if they are antiparallel; here Xq is the coordinate of the interface between the two magnetic layers. At the interface 
between magnetic and nonmagnetic layers (FM/NM) only the longitudinal source term exists, irrespective of the 
alignment of neighboring magnetic layers; it is 



2D 



MiS(x ). 



(A6) 



When two identical magnetic layers are noncollinear there is a transverse source term given by Eq. (A4) as well as a 
longitudinal one Eq. (A3). 

For the multilayered structure depicted in Fig. [l], which models the case studied up till now fo r cu rrent induced 
switching, no two magnetic layers are adjacent so the sole source term that exists is give n by Eq. (A3). In this case 
the boundary condition at the interfaces between adjacent FM/NM layers is given by Eq. (A6). However, as we make 
the assumption that the thickness of the nonmagnetic spacer layer between the two magnetic layers is much smaller 



than we can replace the two sets of FM/NM boundaries by one and use the conditions Eqs. (|A3|) and (A4) 



APPENDIX B 



In this Appendix, we derive the boundary conditions at the interfaces between the layers in Fig. [j]. To achieve 
this goal, we consider a sub-system shown in Fig. ^ which consists of a semi-infinite FM layer with x < with 
the local magnetization M = cos 8e z — sin 8e y , a diffuse interfacial layer I between < x < d 1 with the same local 
magnetization as in the FM layer, and a semi- infinite NM layer for x > d 1 . When d 1 is infinitesimally small, this 
sub-system represents the three FM-NM interfaces in Fig. [I], i.e., between the thick FM and spacer layers, between the 
spacer and thin FM layers (when spatially inverted), and between the thin FM and back NM layers. We assume that 
both spin accumulation and current are continuous at the FM-I and I-NM interfaces, and derive the relation between 
spin accumulation and current at x = with the same quantities at x — d 1 as the thickness of the interfacial layer 
d 1 goes to zero. In this limit the parameters of the interfacial layer, such as \ I m f pl t^, J 1 , Aj, and, most important, 
its resistance ARj remain constant; the latter condition implies that the diffusion constant of the interfacial layer 
~ d 1 as d 1 0. 

We solve Eqs. (|), (|) for the spin accumulation and use Eq. (||) to find spin current in the ferromagnetic and 
interfacial layers. By adopting a set of local coordinates (x, y, z) such that the local magnetization is Mj = the 
spin accumulation and current in the FM layer take the form: 

' mf = 2Re (g 2 exp(^)) 

< m F y = 2/m(G 2 exp(^)) (Bl) 
= Giexp( I ^), 
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FIG. 6: Structure of the interface between the layers at Fig.l. FM is a semi-infinite ferromagnetic layer with the local 



magnetization NL d — cos8e z — sm9e y , I is a diffuse interfacial layer with the same local magnetization as in FM layer, and 
NM is a semi-infinite nonmagnetic layer. 



and 



i F - 

■I III ..r 



i F - 

Jm.y 



i F - 

Jm.z 



-4D F Re^exp(fk)) 
-4^/m(^exp(^)) 

/&- ag °>' g,) Giexp(sH. 



In the interfacial layer, spin accumulation and spin current take the form: 

' ml = 2Re (G 5 exp(-$-j) + 2Re fG 6 exp(-^-) 



ml = 27m (G 5 exp(f )) + 2/m (G 6 exp(-^) 
m\ = G 3 exp( 1 4^) + G 4 exp(-^-), 



and 



3m,x 



Re (|> exp(^)) - Re exp(-^ 
j 1 ^ = -4D^ \lm (ff exp(^)) - Im exp(-^) 



2D I a (l- 11 ') 



G 3 exp( I ^) - G 4 exp(- I ^) 



(B2) 



(B3) 



(B4) 



Here /3, are spin-polarization parameters for the conductivity defined in Sec. 0; 7, 7' are similar parameters for 
the diffusion constant defined as Cj — ^Cl'M.d in the bulk of the ferromagnetic layer, and Dj = ^'D^lsAd for the 
interfacial layer, and = 



1 



(A-)* (a j ,J ) 2 ^ whenA?'«A fl/ 



F.I 



The boundary conditions for the continuity of the spin accumulation and current at the interface between ferro- 
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magnetic and interfacial layer x — take the form: 

2ReG 2 = 2ReG 5 + 2ReG e 
2ImG 2 = 2ImG b + 2ImG e ) 
G\ = G3 + G4, 



(B5) 



and 



'+ + 
-AD$Im{%L) = -4Z>£lm(%^) 

Pj. ~ 2D '^ ) G 1 = Tie - ^%^(G 3 - G 4 ). 



(B6) 



To relate m f (0) to m / ((i / ), and j„(0) to we use the assumption that as the thickness of the interfacial layer 

goes to zero, other parameters of the interfacial layer, such as Xl dl , J 1 , and Xj remain constant, but the diffusion 
constant Dq goes to zero with the same rate as d 1 , so that d 1 /Dq = const. Then, for example, for small d 1 <§; Aj the 
x-component of the spin-accumulation at x — d 1 may be written as 



miid 1 -»■ 0) « 2Re(G 5 + G 6 ) + 2Re ( (G 5 - G 6 )^- 



(B7) 



By comparing this expression with Eqs. (B5) and ( |B6| ), we obtain a relation between the x-components of spin 
accumulation and current at x = and x = d 1 : 



miid 1 -0) = m£ (0)-j£ =(0) 



2D<> 



(B8) 



and similarly, 



mUd 1 



0) = mf (0) - £ ,(0) 



2^' 



(B9) 



m£(d' - 0) = mf(0)+j ' " -j£ (0) — - " , 

2(1 - 77') D Q 2(1 - 77') D J 



(B10) 



In a manner similar to Eq. (B7) the x-component of spin current at x = d 1 may be written as 



G5 — Ge 



2i?e (i(G 5 + G 6 )) 



By comparing this expression with Eqs. (|B5| ) and (B6), we find the continuity condition for the a;-component of spin 
current: 



3Lx(d I ->0)=£A0)-rn*(0) 



and, similarly, 



(Bll) 



jLyid 1 - 0) = j£ iff (0) + m£(0) 



(B12) 



and 



im,6ar,(^ - 0) = j£ )Z -(0) - mf (0) — . 



(B13) 
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With these relations we can now obtain the boundary conditions at the three interfaces in the multilayered structure 
depicted in Fig. |l]. 

By using the conditions Eqs. (B8)-(B13), the boundary conditions at the interface between the thin (first) ferro- 
magnetic and non-magnetic (N) layers of the structure shown in Fig. ^ at x — t F may be written immediately, since 
in the thin FM layer the local coordinate system (x, y, z) coincides with the global axes [x, y, z); we find 

' mg(t F ) - m x 1] {t F ) = -rj£l(tF) 



m y (*f) - m ( y\t F ) = -rj ( „l]y(t F ) 
[ mf (t F ) - m£\t F ) = rj eJ ^r ~ rj&fo)^ 



(B14) 



and 





jm y x (t F ) 


= -m y l \t F )£f 




- jm]y(t F ) 






-3 { rk ] At F ) 


= -m z 1 \t F )£-, 



(B15) 



where r = d 1 /2Dq. Similarly, the boundary conditions at the interface between the non-magnetic spacer (S) and the 
thin (first) FM layer at x = take the form: 



m^(0)-m^ 1) (0)=ni 1 , ) 2/ (0) 
^ mf (0) - mi X) (0) = -rj eJ 



(B16) 



77' 



and 



f jl*(0)-j£k0) = m«(G) 



d'.T 



(B17) 



(!)/ 



Note that spin-current conservation condition at the interfaces, which means that there are no torques acting at the 
interfaces, is due to the infinitely small thickness of the interfacial layers d 1 — ► 0. To write the boundary conditions 
at the interface between the thick FM and NM spacer layers at x = 0, we have to change from the local coordinate 
system (x, y, z), related to the magnetization direction in the thick FM layer, to the global (x,y,z) system. Any 
vector a will be transformed according to the following rule: 



ay = a y cos 8 + a z sin 8 
a z = —a y sin 8 + a z cos ( 



(B18) 



By applying this transformation to the conditions, Eqs. (B8)-( B13 ), we obtain the following boundary conditions at 
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the interface between the thick (second) FM and non-magnetic spacer (S) layers: 

mf(0)-m( 2) (0) = ~rji 2) ,(0) 



and 



rriy(0) - my(0) = -rj 



I-77 



7 sin ( 



(2) /-m I-77' cos 2 



1-77' 



+rj$ B (0) sin( 



(0) - mi 2) (0) = rj e ^ cos 9 + r (0) sin 9 cos 0^ 



77 
1-77' 

(2) 



-rj&(0)*=? 



1 — 77 / sin 2 



-77' 



(B19) 



i£,*(0) " ji 2 k0) - -m< 2 >(0) cos^ - mi 2) (0) sin^ 
jLy (0) - „ (0) = mi 2) (0) cos ^ - m< 2) (0) sin 2 9 ^ 

+rni 2) (0)sin6'cose'^- ( B 20) 
Jm,*(0) " Jmko) = mi 2) (O)sin0^ + m 2y (0) sin^cos^^- 
-mi 2) (0) cos 2 #4^, 

Note that as the thickness of the interfacial layer d J go es to z er o, for diffuse scattering we considered one produces 
large discontinuities in the spin-accumula tion (E qs. (B14), ( | B16 ), ( B19| )) proportional to finite r = d 1 /2Dq, but small 
discontinuities in the spin-currents (Eqs. ( B15| ), ( B17 ), QB20[ )) proportional to d 1 , because J/ does not increase and r/j 
does not decrease as d 1 —> 0. In our picture finite thickness of the interfacial layer is essential for torque production 
at the interface. As we consider infinitely small interfacial thicknesses d 1 — > 0, we obtain spin-current conservation 
conditions at each interface: 



(B21) 



Jm(0) 



(0), 



(B22) 



and 



£ ) (o)=£(o). 



(B23) 



By eliminating m s (0) and j^(0) from Eqs. (B16), ( B19 ), ( B22 ), and (B22), we finally obtain the boundary conditions 
at the interface between thick (second) and thin (first) FM layers at x — (of course there is the NM spacer in-between, 
however its thickness tjy is irrelevant for these boundary conditions as long as tjy -C A^ ; ): 



m£ x) (0) - mi 2) (0) = -2r^(Q) 
TO™(0)-TO^ 2) (0) = -ri er ^ 



-77 



7 smf 



.(1) |- n ^ 2-7 7 '(l+cos 2 9) 
1 Jm,V\ u ) I-77' 



+rjm]z (0) sin 9 cos 6*y 



77 



(B24) 



mi 1} (0) - to^ j (0) = rj eT7 ^ ? (l + cos 9) + rj^(0) si 

_„,•(!) 2-77' sin 2 9 
ijm,z\y) 1-77' ' 



,( 2 



1 77 
1-77' 
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and 



£ ) (0)=j^ ) (0) 



(B25) 



Finally, we show that parameter r = d 1 j2D\ is proportional to the interface resistance ARi found from CPP 
transport measurements tL!3. By considering the expression (Q) for the electrical current in the interfacial layer, 
and the assumptions that Dq ~ d 1 and \l dl remains constant as the thickness of the interfacial layer d 1 — > 0, we 



find AR,! = d I /2C^, or r 



2D, 



ARijS-. The parameters Gq and Dq may be related via Einstein's relation 



Cj = e 2 Nj(eF)Di , so that the parameter r takes the form: 

r = AR I e 2 Nl(e F ) 



1 - 7 " 2 
1 — 77' ' 



(B26) 



where e is the electron charge, iVg (ep) is the density of states at the interface at Fermi energy, and 7" is the spin 
polarization parameters for the density of states at the interfaces which is defined as N/ = r )" Nq M<j. 



APPENDIX C 



We solve the Eqs. @), (§) for the spin accumulation in each of three layers, and find spin currents using Eq. (j^). 
In the thick ferromagnetic layer, spin accumulation and current take the form 



,(2) 



2Re fG 2 exp(^) 



,( 2 ) - 



rn 



(2) 



27m ^G 2 exp(j^)J cos# — G\ exp(p^) sin# 
27m (g 2 exp(f-j) sind + Gi exp(^) cos0, 



and 



;(2) _ 



—ADnRe 



(t ex p(ir) 



& = -0je sin - 4D Im exp(^)) cos 8 
-Gi exp(-rf-) sin 8 



2Ai(l-/3/3') , 



t = /3j e cos0 - 4D Im exp(^)) sir 



2D (l-/3/3') 



Gi exp(-j-f-) cos( 



In the thin ferromagnetic layer 



2Re (G 5 exp(-^)) + 2Re (g 6 exp(^) 
2Jm (g 5 exp(-^-)) + 21m (g 6 exp(2^)) 
G 3 exp(-^)+G 4 exp(^), 



and 



Re (Gl exp(-^)) - i?e exp(^)) 
j m (G5 ex p(-^)) - 7m exp(^) 



Jm,z — 



G 3 exp(~ I ^)-G 4 exp(^) 



(CI) 



(C2) 



(C3) 



(C4) 
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where 1+ 1 



and 



, and \ f dl is spin-diffusion length in FM layer. In the non- magnetic layer, 

x — t F > 



m = A exp(- 



A sdZ 



-)• 



2D» 



(C5) 



(C6) 



To obtain the 12 unknown constants A x , A y , A Zl G*i, ReG i^ ImG i-, G3, G4, ReG§, ImG$, ReG§, ImGe, we use 
the boundary conditions (see Appendix || Eqs. ([B14| ), ( |B24| ), ( |B2l| ), and (|B25|)): 



and 



m f (*f) - ml (*f) = -rim,»(tf) 
m f - m { y\t F ) = -rj$ y (t F ) 
[ m?(t F ) - m?\t F ) = rj eT ^r - rj^tp)^, 

m£ x) (0) - mi 2) (0) = -2rj£l(0) 
m«(0) - (0) = -rie^ sinfl - ^(0)^12^^ 
+rjl 1 , ) 2 (0) sin flcose^ 

,( 2 



m^(0) - m^(0) = rje^^l + cosfl) + rj^(O) si 

..•P) fn'l 2 -77' sin 2 9 



1-77' 



(C7) 



(C8) 



ti(t F )=^>(t F ) 



(C9) 



£ ) (0)=£ ) (0) 



(CIO) 



where the parameter r is proportional to the interface resistance ARi , r = ARie 2 No(l — 7" 2 )/(l — 77') , e is the 
electron charge, Ao is the density of states at the interface, 7, 7', 7" are the spin polarization parameters for the 
conductivity, diffusion constant, and density of states at the interfaces (see Appendix [bJ). The other six boundary 
conditions come from the conservation of spin current at the interfaces. 
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